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Abstract A linear elastic fracture mechanics model for cal¬ 
culation of the splitting strength of dowel-type fastener 
joints loaded perpendicular to grain (Van der Put/Leijten 
model) has previously been presented, and now forms the 
basis for design in Eurocode 5. The original Van der Put/ 
Lcijten model was derived using a number of simplifying 
assumptions, e.g., that the normal forces in the cracked 
parts of the beam can be ignored, leading to a solution that 
does not involve the effect of an initial crack. In the present 
article an extended version of the Van der Put/Leijten 
model is derived without any simplifying assumptions, and 
it is shown that the original Van der Put/Leijten model 
appears as a special case, namely by assuming that only 
contributions from shear deformations are significant. The 
model presented here involves the effect of an initial crack 
and may be characterized as a generalized linear elastic 
fracture mechanics model. Results of tests showing the in¬ 
fluence of initial cracks of various lengths are presented and 
compared with the predictions. 

Key words Dowel joints ■ Perpendicular to grain ■ Split¬ 
ting ■ Fracture mechanics ■ Crack length 


Introduction 

A dowel-type fastener joint loading a beam perpendicular to 
grain may fail in a ductile manner, characterized by 
bending of the fastener and/or embedment of the fastener 
into the wood, or cause brittle failure characterized by split¬ 
ting of the beam. The ductile failure modes are well under¬ 
stood and can be accurately predicted by the European yield 
model 1 (or extended theories), which now forms the basis of 
design of dowel-type fastener joints in major design codes. 
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For brittle failure modes, no simple theory suitable 
for implementation in design codes has yet gained wide 
acceptance. Recently, however, a number of simple theo¬ 
retical models based on fracture mechanics have been pro- 
posed. 21 A large number of test data compiled from the 
literature was analyzed, 4 and the models seem to be able to 
predict the splitting failure fairly well (at least if introducing 
appropriate empirically determined effectiveness factors). 
The Van der Put/Leijten model 2-4 now forms the basis of 
design in Eurocode 5. 7 

The original Van der Put/Leijten model was derived 
based on a number of simplifying assumptions, among 
others, the disregard of normal forces in the cracked part of 
the beam, leading to a very simple formula readily appli¬ 
cable to practical design. An extended version of the Van 
der Put/Leijten model was presented, 8 which takes into ac¬ 
count the normal forces in the cracked parts of the beam, 
and does not resort to simplifying assumptions. The result¬ 
ing formula for the splitting strength remains fairly simple, 
and the original Van der Put/Leijten solution appears as a 
special case if only the contribution from shear deforma¬ 
tions is considered or if the crack length is assumed to be 
zero. 

A basic limitation for application of linear elastic frac¬ 
ture mechanics (LEFM) is that the size of the fracture pro¬ 
cess zone must be small when compared with the length of 
the crack. This condition is not fulfilled in the Van der 
Put/Leijten model (at least for short cracks). A simple and 
widespread modification of LEFM to overcome this prob¬ 
lem is to increase the assumed crack length. 9 Such a modifi¬ 
cation is easily accomplished using the model derived in the 
present article. 


Strength analysis 

Linear clastic fracture mechanics 

For a linear elastic body loaded by a single load, P, the crack 
propagation energy release rate, <g , is given by 10 







481 


P * 2 dC 
2i> dfl 


( 1 ) 


where £> is the width of the body, a is crack length, and C is 
the compliance given by 


C = - (2) 

P 

where b is the deflection of the loading point. 

A crack starts propagating when the energy release rate 
reaches the critical value, (g z . Assuming static or quasi static 
conditions and no energy dissipation outside the fracture 
region, (g z is equal to the fracture energy, Gj f , of the 
material. 

The crack propagation load, P c , of the body is thus 
given by 


P = 



( 3 ) 


A simply supported beam loaded perpendicular to grain 
by a single joint is shown in Fig. 1. It is here assumed 
that the fastener is sufficiently stiff to ensure that the 
crack propagates along the grain simultaneously through 
the entire width of the beam. 

The cracked beam is modeled using ordinary beam 
theory. The static model used is shown in Fig. 2 for a 
symmetrical crack of length a = 2ph, h being the 
beam depth. The span of the beam is 2 L, and the load 



( i- a yi 

ah 


Fig. 1. Geometry, load, and support conditions 


| L-ph j ph ph | Lrflh j 



Fig. 2. Static model for cracked beam 


applied at midspan is denoted P. The edge distance is 
denoted ah. 

Beam theory should formally not be expected to render 
useful solutions for short cracks (e.g., 2 fth < ah), but is 
attempted here anyway. The beam structure shown in Fig. 2 
is three times statically indeterminate. The unknown vari¬ 
ables N 2 , V 2 , and M-, are chosen as the sectional forces at the 
center of beam 2, and are found to be 



The deflection at the loading point, b , is given by 


b = P 


-—I L + /3h 
5 GA' 1 


1 — a 




3 L(L - Ph) + (Ph) 2 - -(2 L - phf{ 1 - a 3 ) 


( 5 ) 


where E is modulus of elasticity (MOE), G is shear 
modulus, A is cross-sectional area, and I is moment of 
inertia of the uncracked beam cross section (width b, depth 
h). 

By means of Eqs. 3 and 5, P c is found to be 



ahG(g t 


-“>+If (!)<-«■) 


( 6 ) 


ft is noted that P c is a decreasing function of p, i.e., stable 
crack growth is not possible, and the crack propagation load 
according to the present model is thus equal to the failure 

load causing catastrophic failure. 

Zero crack length (P —» 0) or only considering the contri¬ 
bution to the deflection from the shear force (E —» °o or 
simply disregarding the second term in Eq. 5) leads to the 
failure load, / J c0 . which is the original Van der Put/Leijten 
solution 2 4 


P c0 = 2bJjG&h-?~ (7) 

V 3 1 — a 

The extended Van der Put/Leijten solution as given by Eq. 
6 may then also be expressed as 


5 G 1 - a 3 * * * f pY 

2 E 1 — a\aj 

A general LEFM solution for a beam with a crack, based 

on simple engineering beam theory, has previously been 

presented. 11 The left side of the crack (Figs. 1 and 2) is 

shown in Fig. 3. The sectional forces at the crack tip are 

given by Eq. 4 and by 
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Fig. 3. Sectional forces at crack tip 


Nr = -N 2 , V 3 = 0 

V 3 = -P, V 3 = -P 
2 2 

M, = i P{L - Pti) + |WV 2 - M 2 , 

M 3 = |P(L - Ph) 

The crack propagation load, P c , is given by 11 


( 9 ) 
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where A = bh , A 2 = aA, A 2 = (1 — a)A , I = bb?l 12, / 3 = a 3 /, 
/ 2 =(l-a) 3 /. 

Inserting the sectional forces as given by Eqs. 4 and 9 in 
Eq. 10, P c is found to be given by Eq. 6. 


Generalized linear elastic fracture mechanics 


A basic limitation for application of LEFM is that the size of 
the fracture process zone must be small as compared with 
the length of the crack. A simple modification of LEFM to 
overcome this problem is addition of a crack length, Aph. 
An expression for estimation of such an additional crack 
length for mode I fracture of wood based on nonlinear 
fracture mechanics calculations using the finite element 
method was given by 9 


Aph 


0.2 ^ f V£G 

J 2 


( 11 ) 


where / t is the tensile strength perpendicular to grain. 

The failure load of a structure is, according to Eq. 8, 
obtained for a minimum value of the crack length. The 
minimum crack length to be inserted in Eq. 8 may be esti¬ 
mated as ph = Aph + dl 2, d being the diameter of the dowel 
hole, leading to 



/. [MPa] 


Fig. 4. Influence of perpendicular-to-grain tensile strength 


Pc = «/ J c0, 

L , 5 1-a 3 G 2 fo.2g t | d V (12) 

] 2 1 - a (ah ) 2 1 f t 2 2 JeG J 

According to Eq. 11, the tensile strength perpendicular 
to grain has a significant influence on the additional crack 
length. Figure 4 shows an example of the influence of the 
tensile strength on the failure load using h = 300 mm, ah = 
80 mm, d = 20 mm, E = 12 700 MPa, G = 870 MPa, § f = 
0.25 N/mm [material properties for laminated veneer lum¬ 
ber (LVL) 5 ]. 

The perpendicular-to-grain tensile strength, which is sel¬ 
dom used and reported in the literature, is highly volume 
dependent, and an appropriate value for use in the present 
model is not obvious. 

The splitting strength of LVL beams loaded perpendicu¬ 
lar to grain by dowel joints has previously been tested. 6 The 
specimens used for determination of the perpendicular-to- 
grain tensile strength for LVL were 40 X 70 X 280 mm 3 
(70mm in parallel-to-grain direction, 280 mm perpendicular 
to the veneer plane), resulting in a reported tensile strength, 
f = 0.89 MPa. 3 The size of the tensile test specimens may be 
suspected of being somewhat too large for the present pur¬ 
pose, thus leading to a somewhat low tensile strength. 

For glulam and solid timber, the original Van der Put/ 
Leijten solution as given by Eq. 7 may be sufficiently accu¬ 
rate for practical use, but for materials with very low per- 
pendicular-to-grain tensile strength it may significantly 
overestimate the failure load. 


Experimental 

The Van der Put/Leijten model (Eq. 7) has previously been 
compared with a considerable number of tests on glulam 
beams. 412 The tests presented in this article were conducted 
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primarily to evaluate the appropriateness of the modifica¬ 
tion factor, u, as given by Eq. 8. 

Methods and materials 


Results and discussion 

A so-called fracture parameter, C l5 was proposed and de¬ 
rived from plate-joint specimen tests 13 


Test series 1 

Glulam beams without finger joints were tested in three- 
point bending as shown in Fig. 1. One 14-mm dowel was 
used in a 15-mm hole. All beams had a span of 1200mm, 
width b = 25 mm, and edge distance h e = 40 mm. Two 
different beam depths, h = 100 mm (a = 0.4) and h = 
200mm (a = 0.2), were tested. 

Initial crack lengths /3h = 0, 20, 40, 80, and 160mm (i.e., 
/3 la = (3h lh e = 0,0.5,1,2, and 4) were tested. Six specimens 
were tested for each condition. 

The glulam was made of Japanese cedar ( Cryptomeria 
japonica), lamella thickness 30 mm, and all lamellae were 
the same grade (machine graded). MOE was measured on 
the five glulam beams (100 X 200 mm 2 cross section, length 
3000 mm) from which the specimens were cut by measuring 
the longitudinal vibration frequency, resulting in E = 
7530 MPa. Moisture content (MC) at the time of testing was 
10.5% and density at the given MC was 372 kg/m 3 . 

Furthermore, tests were conducted on plate-joint speci¬ 
mens (nine specimens were tested) as shown in Fig. 5 in 
order to derive the necessary material properties. 

Test series 2 

Test series 2 was conducted to evaluate the influence of the 
edge distance, h e , of the plate joint on the fracture param¬ 
eter as given by Eq. 13 (or the apparent fracture energy 
derived here from). 

A total of 40 plate-joint specimens (dimensions as shown 
in Fig. 5) were cut from five glulam beams, which had a 
mean MOE of 5670MPa. Three different edge distances 
were tested: /? e = 20mm (10 specimens), /7 e = 40 mm (20 
specimens), and h e = 60 mm (10 specimens). 



(13) 


where P c is the recorded failure load. Equation 13 appears 
from Eq. 7 by assuming h —» °o. 


Test series 1 

The plate-joint tests resulted in C l = 9.38 ± 1.40 N/'mm 2 , or 
if calculating the fracture energy using G = £718: <g t = 0.129 
± 0.039 N/mm. 

Figures 6 and 7 show the modification factor, ft = PJP C0 , 
as obtained from the beam tests and the theoretical solution 



/3/a 

Fig. 6. Modification factor, fi, as function of normalized crack length 
for h = 100 mm. Solid line , Eq. 8; filled circles , knots in vicinity of both 
crack tips 



Units: mm 



Fig. 7. Modification factor, fi, as function of normalized crack length 
for h = 200 mm. Solid line, Eq. 8; filled circles , knots in vicinity of both 
crack tips 


Fig. 5. Plate-joint specimen 
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P la 

Fig. 8. Failure load as function of normalized crack length for h = 
100 mm. Solid line, Eqs. 7 and 8; filled circles, tested mean value. Bars 
indicate standard deviation 



P/a 

Fig. 9. Failure load as function of normalized crack length for h = 
200 mm 


(Eq. 8). For the test results, the mean value of the failure 
load obtained for /3h = 0 (or in reality /3h = 7.5 mm due to 
the dowel hole) has been used for P c0 . (Note: for 
h = 100 mm as well as for h = 200 mm, one excessively high 
failure load has been discarded in calculation of Pco)- 
Note that the only material property involved in « is the 
ratio GIE. 

In Figs. 8 and 9, the tested and calculated (Eqs. 7 and 8) 
failure loads are shown. Figures 6 and 7 indicate that the 
modification factor, u, as given by Eq. 8 predicts an influ¬ 
ence of the crack length in good agreement with test results. 
Figures 8 and 9 indicate that the original Van der Put/ 
Lcijten formula as given by Eq. 7 is in good agreement with 
test results if the so-called fracture parameter 13 (or apparent 
fracture energy) is determined by means of plate-joint 
specimens (Fig. 5). 

According to Eq. 6, stable crack growth is not possible; 
catastrophic failure is deemed to occur once the crack starts 
propagating. This prediction is, however, not in accordance 


Table 1. Fracture properties determined by plate-joint test 


h c (mm) 

Cj (N/mm 3 ' 2 ) 

Cjt (N/mm) 

20 

9.68 ± 1.29 

0.181 ± 0.048 

40 

9.95 ± 1.56 

0.193 ± 0.062 

60 

9.31 ± 1.85 

0.171 ± 0.071 

Mean ± sd 

9.72 ± 1.56 

0.185 ± 0.060 


Q t based on G = El 18 


with the experiments, in which cases of crack growth were 
observed during increasing load. The loads given in the 
figures are the catastrophic failure loads; no attempts were 
made to measure the crack lengths at failure. 

Test series 2 

It may be suspected that the good agreement between pre¬ 
dicted and tested failure loads as shown in Figs. 8 and 9 
obtained by use of the plate-joint specimens is due to the 
fact that the same constant edge distance ( h e = 40 mm) was 
used in all plate-joint specimens and beam specimens. Test 
series 2 was conducted in order to evaluate the influence of 
h,, on the fracture parameter, Cj. 

The plate-joint specimens resulted in the values of the 
fracture parameter, Cj, and corresponding apparent frac¬ 
ture energy, <g { , as given in Table 1. Table 1 shows that the 
fracture properties derived from plate-joint tests using Eq. 
13 are constant in the tested range of edge distances. 

It may be worthy of notice that the fracture parameter, 
C\, obtained for the two test series (E = 7530 MPa and E = 
5670 MPa) only varies insignificantly. However, if G is as¬ 
sumed to be a fixed ratio of E, Eq. 13 leads to significantly 
different values of the fracture energy for different E. Thus, 
if as usual a constant value of the fracture energy is used 
regardless of the MOE, the model presented here (and the 
Van der Put/Leijten model) may give somewhat poor pre¬ 
dictions for some values of the MOE. 

Such large variation in fracture energy due to variations 
in MOE as reported here is not in accordance with experi¬ 
ences using double cantilever beam (DCB) specimens for 
determination of the fracture energy. It is unclear at present 
whether the observed difference in fracture energy is due to 
an improper estimation of G as a constant fraction of E, or 
due to improper model assumptions, although the latter is 
the most likely. At any rate, use of the fracture parameter, 
Cj determined from plate-joint tests using Eq. 13 seems to 
make theoretical predictions as given by Eq. 6 in good 
agreement with test results. 

Although it would be more satisfactory to be able to use 
common values of shear modulus and fracture energy 
rather than resorting to introducing a new property like C u 
which needs to be determined by some special tests, 
this procedure is by no means unique. For instance in the 
treatment of glued lap joints, 14 glued-in steel rods, 13 and 
glued-in hardwood dowels, 16 it is commonly accepted 
that special tests are needed for determination of the bond 
line parameters, the values of which may be orders of 
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magnitude different from those determined directly on the 
adhesive. 

Testing standards are available and commonly used for 
testing of material properties (embedding strength and 
foundation moduli) needed as input in calculation models 
for prediction of ductile failure modes in dowel-type fas¬ 
tener connections. However, no testing standards seem 
available for testing of properties relevant to brittle splitting 
failure modes, ft seems reasonable to suggest such standard 
tests, and because the plate-joint specimens are very easy to 
make and the tests very easy to conduct, this specimen is 
recommended as a standard test. 


Conclusions 

A linear elastic fracture mechanics model for calculation of 
the splitting failure load of beams loaded perpendicular to 
grain by dowel-type fastener connections was presented. 
The new model is based on an approach similar to the 
previous Van der Put/Leijten model, 2 but includes normal 
forces in the cracked parts of the beam. The Van der Put/ 
Leijten solution appears as a special case, namely for zero 
crack length or by only considering contributions from 
shear deformations. 

Like the Van der Put/Leijten model, the presented 
model is two-dimensional, i.e., it assumes that the fasteners 
are sufficiently stiff to ensure that crack propagation along 
the grain occurs simultaneously through the entire width of 
the beam. Especially for joints with large relative edge dis¬ 
tances (a values), this condition may not always be fulfilled 
because the dowels are designed sufficiently slender so as to 
enable ductile bending at the ultimate limit state. 

Because the new model includes the effect of crack 
length, it is easily extended from LEFM to generalized lin¬ 
ear elastic fracture mechanics by adding the finite size of the 
fracture process region to the crack length. Tests were con¬ 
ducted to evaluate the influence of crack length, and good 
agreement was found between theoretical predictions and 
test results. Furthermore, tests were conducted to derive 
appropriate material properties needed as input param¬ 
eters. It is strongly recommended that the so-called plate- 
joint specimen be used as a standard test specimen for 
determination of relevant fracture properties. 


References 

1. Johansen KW (1941) Fors0g med trarforbindelser. 
Bygningsstatiske Meddelelser, argang XII 


2. Van der Put TACM, Leijten AJM (2000) Evaluation of perpen¬ 
dicular to grain failure of beams caused by concentrated loads of 
joints. In: Proceedings of the International Council for Research 
and Innovation in Building and Construction, CIB-W18 Meeting 
Thirty-three, Delft, The Netherlands, paper no. 33-7-7 

3. Leijten AJM, Jorissen AJM (2001) Splitting strength of beam 
loaded by connections perpendicular to grain, model validation. In: 
Proceedings of the International Council for Research and 
Innovation in Building and Construction, CIB-W18 Meeting 
Thirty-four, Venice, Italy, paper no. 34-7-1 

4. Leijten AJM (2002) Splitting strength of beams loaded by connec¬ 
tions, model comparison. In: Proceedings of the International 
Council for Research and Innovation in Building and Construc¬ 
tion, CIB-W18 Meeting Thirty-five, Kyoto, Japan, paper no. 
35-7-7 

5. Gustafsson PJ, Larsen HJ (2001) Dowel joints loaded perpendicu¬ 
lar to grain. In: Aicher S, Reinhardt HW (eds) Joints in timber 
structures. Proceedings of the International RILEM Symposium, 
Stuttgart, Germany, pp 577-586 

6. Jensen JL, Gustafsson PJ, Larsen HJ (2003) A tensile fracture 
model for joints with rods or dowels loaded perpendicular-to- 
grain. In: Proceedings of the International Council for Research 
and Innovation in Building and Construction, CIB-W18 Meeting 
Thirty-six, Colorado, USA, paper no. 36-7-9 

7. EN 1995-1-1 Eurocode 5 (2002) Design of timber structures. Part 
1-1: General rules and rules for buildings. Comite Europeen de 
Normalisation, Brussels, Belgium 

8. Jensen JL (2003) Splitting strength of beams loaded by connec¬ 
tions. In: Proceedings of the International Council for Research 
and Innovation in Building and Construction, CIB-W18 Meeting 
Thirty-six, Colorado, USA, paper no. 36-7-8 

9. Gustafsson PJ (2003) Fracture perpendicular to grain - structural 
applications. In: Thelandersson S, Larsen HJ (eds) Timber engi¬ 
neering. Wiley, Chichester, pp 103-130 

10. Hellan K (1984) Introduction to fracture mechanics. McGraw-Hill, 
New York, pp 49-72 

11. Petersson H (2002) Energy release rate analysis. In: Gustafsson PJ 
(ed) Fracture mechanics models for strength analysis of timber 
beams with a hole or a notch - a report of RILEM TC-133. Report 
TVSM-7134, Lund University, Sweden 

12. Ballerini M, Giovanella A (2003) Beams transversally loaded by 
dowel-type joints: influence on splitting strength of beam thickness 
and dowel size. In: Proceedings of the International Council for 
Research and Innovation in Building and Construction, CIB-W18 
Meeting Thirty-six, Colorado, USA, paper no. 36-7-7 

13. Yasumura M (2002) Determination of fracture parameter for 
dowel-type joints loaded perpendicular to wooden grain and its 
application. In: Proceedings of the International Council for 
Research and Innovation in Building and Construction, CIB-W18 
Meeting Thirty-five, Kyoto, Japan, paper no. 35-7-9 

14. Gustafsson PJ (1987) Analysis of generalized Volkersen-joints in 
terms of non-linear fracture mechanics. In: Verchery G, Cardon 
AH (eds) Mechanical behaviour of adhesive joints. Pluralis, Paris, 
pp 323-338 

15. Gustafsson PJ, Serrano E (2001) Glued-in rods for timber struc¬ 
tures - development of a calculation model. Report TVSM- 
3056, Division of Structural Mechanics, Lund University, Lund, 
Sweden 

16. Jensen JL, Koizumi A, Sasaki T, Tamura Y, Iijima Y (2001) 
Axially loaded glued-in hardwood dowels. Wood Sci Technol 
35:73-83 




